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We present two recently obtained solutions of the Einstein equations with spherical 
symmetry and one additional Killing vector, describing colliding null dust streams. 

1 Introduction 

Null dust is a model for the incoherent superposition of waves with random 
polarizations and phases, moving in a single direction with the speed of light. 
It also describes the high frequency (geometrical optics) approximation for 
any type (including gravitational) radiation. Null dust is characterized by the 
energy-momentum tensor T ab — pl a l b . Neither the energy density p nor the 
null vector I are uniquely determined and the former can be absorbed in the 
latter by an appropriate rescaling. Various exact-, solutions of the Einstein 
equations in the presence of null dust are knowuJil&tj, but certainly the most 
famous of them is Vaidya's shining star solutionQ. 

The energy-momentum tensor describing the collision of two streams of 
null dust (with propagation vectors k and I pointing in opposite spatial direc- 
tions) is a sum of the two single null dust energy-momentum terms: 

T ab = ri b + k a k b ^ 

The Einstein equations for the matter source (Q) in the plane symmetric case 
were solved by Letelier a. The cylindrically symmetic case was studied by 
Letelier and WangH. Progress in the numerical treatment of the spherically 
symmetric case was achievecLby Holvorcem, Letelier and WangQ. Under the 
assumption of staticity DateQ integrated half of the relevant equations. The 
solution was found by Kramer □ and by usE3, ours containing two additional 
parameters. Later, assuming homogeneity in addition to spherical symmetry, 
we have obtained a Kantowski-Sachs type solution: a closed Universe filled 
with colliding streams of radiation M. The interpretation of such spherically 
symmetric solutions as regions of space-time containing both escaping and 
backscattered gravitational radiation was proposed by Poisson and Israel 113. 



Recently Kramer obtained static and stationary solutions describing the colli- 
sion of cylindrical null dust beamsli3. 

Here we present briefly the spherically symmetric static and homogeneous 
solutions. 



2 The spherically symmetric static and homogeneous solutions 

The detailed derivation of these solutions was presented elsewhere 00. The 
line element including both the homogeneous case (c = 1) and the static case 
(c = —1) can be given in the concise form: 

ds 2 = — — (dZ 2 - 2R 2 dL 2 ) + R 2 dn 2 . (2) 
cCR 

Here R(L) is a transcendental function (depending on the error function): 

2 f L 2 

cCR = e L2 -2L<f> B (L), $ B (L) = B + dx > , (3) 

while B and C > are constants of integration. 

One may think of Z as being a time coordinate t in the static case and a 
radial coordinate r in the homogeneous case. Conversely, the coordinate L has 
the meaning of time in the homogeneous case and is a radial coordinate in the 
static case. However, for this interpretation to hold, a study of the sign of the 
expression (||) of cCR is necessary, more precisely in both cases_R_> should 
hold. This requirement follows also from the energy conditions li9'Eil. 

First of all let us remark that there is a true singularity at R = 0. This 
can be seen for example from the ill-behaveness at R = of the Kretschmann 
scalar 1Z a bcdR- abcd oc R~ 6 . Thus the domain of the coordinate L is bounded by 
its value(s) Lq at R = 0. The constant B is related to Lq as: 

e L 2 o rlo 



B - X(L ) = — -Je*dx. (4) 

For each value of the constant B there are two values Lq± , a positive and a 
negative one (Fig.l). These are situated symmetrically with respect to the 
origin only for B = 0. 

Second, we express the function cCR in the alternative form 

cCR = 2L \x(L) - x(L )] . (5) 

Because x(L) is monotonically decreasing: d\jdL = —e L2 /2L 2 < 0, the do- 
mains of positivity of R are Lo- < L < io+ m the homogeneous case and 
L < Lq- together with L > L 0+ in the static case. 



2 




-10 ; 



Figure 1: The curves B = x(-^0±)> corresponding to the singularities R = 0, devide the 
horizontal axis L into domains where R > 0, < 0, > (static case) or R < 0, > 0, < 

(homogeneous case). 



3 The interpretation of the solutions 

In the homogeneous case the time coordinate L is bounded both from above and 
below, the boundaries being singularities. The radius function R is increasing 
for L 6 (Lq-,L*), then decreasing for L e (L* , L +), jiihere L* is defined by 
$b(L*) = 0. The solution represents a closed Universelii. From the symmetry 
group it is obvious that this is a Kantowski-Sachs type solution, difficult to 
relate to any physical situation. 

We have analyzed in detail the static solution for L > Lq + . The metric 
is not asymptotically flat and close to the axis R = (which is the singularity) 
its mass decreases to negative values. Thus we run again into interpretational 
difficulties. However, cutting the interior region along a timelike hypersurface, 
we can glue the solution to some physically meaningful matter. We have 
discusscd-the junction with a generic static interiortHI by applying the Darmois- 
IsraeloEj matching procedure. _Ln particular the junction algorithm with the 
interior Schwarzschild solution tj has fixed the free parameters B and C of 
the colliding null dust. The naturaLexteriors are the incoming and outgoing 
Vaidya regions, which were matchedlij along null hypersurfaces to our solution, 
by applying the Barrabes-Israel junction procedurel!3. 

In the range L < Lq- we have just an other copy of the static solution 
with parameters —B and C, and —L as the radial variable. Indeed, the metric 
(^) is invariant under the interchange (L, B) — > (-L, —B). 

These solutions can be reinterpreted as anisotropic fluids with radial pres- 
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sures equal with their energy densities, and no pressures in tangential-directions 
to the spheres. An other approach we already presented elsewhere til is based 
on the connections with two-dimenssional dilatonic models. In this context 
MikovicEa has given a solution in the form of a perturbative series in powers 
of the outgoing part of the energy-momentum tensor. 

Acknowledgments 

This research has been supported by the OTKA grant D23744 and by the 
Hungarian State Eotvos Fellowship. 

References 

1. D. Kramer, H. Stephani, M.A.H. MacCallum, and E. Herlt, Exact Solu- 
tions of Einstein's Field Equations (Cambridge University Press, Cam- 
bridge, England, 1980). 

2. D. Kramer, in Rotating Objects and Relativistic Physics, Proceedings, El 
Escorial, Spain 1992 , edited by J.J.Chinea and L.M. Gonzalez-Romero 
(Springer Verlag, Berlin, 1993). 

3. J. Bicak and K. Kuchaf, Phys. Rev. D 56, 4878 (1997). 

4. P.C. Vaidya, Proc.-Indian Acad. Sci, Sect. A 33, 264 (1951). 

5. P.S. Letelier, Phys. Rev. D 22, 807 (1980). 

6. P.S. Letelier and A. Wang, Phys. Rev. D 49, 5105 (1994). 

7. P. R. Holvorcem, P.S. Letelier and A. Wang, J. Math. Phys. 36, 3663 
(1995). 

8. G. Date, Gen. Relativ. Gravit. 29, 953 (1997). 

9. D. Kramer, Class. Quantum Grav. 15, L31 (1998). 

10. L. A. Gergely, Phys. Rev. D 58, 084030 (1998). The variables R,r were 
denoted there by r, r* . 

11. L. A. Gergely, submitted. 

12. E. Poisson and W. Israel, Phys. Rev. D 41, 1796 (1990). 

13. D. Kramer, Class. Quantum Grav. 15, L73 (1998). 

14. G. Darmois, in Memorial des Sciences Mathematiques (Gauthier-Villars, 
Paris, 1927), Fascicule 25, Chap. V. 

15. W. Israel, Nuovo Cimento B XLIV, 4349 (1966). 

16. K. Schwarzschild, Sitz. Preuss. Akad. Wiss., 424 (1916). 

17. C. Barrabes and W. Israel, Phys. Rev. D 43, 1129 (1991). 

18. A. Mikovic, Phys. Rev. D 56, 6067 (1997). 



4 



